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The  Turbule  Ensemble  Model  (TEM)  was  developed  to  handle 
acoustical  scattering  from  anisotropic  inhomogeneous  turbulence.  A 
turbule  is  a  localized  atmospheric  inhomogeneity.  TEM,  then, 
represents  a  turbulent  region  by  a  collection  of  turbules  of  different 
sizes.  Since  acoustic  sources  and  sensors  are  omnidirectional,  the 
scattering  volume  of  the  TEM  region  is  ill  defined.  Scattering 
properties  of  individual  turbules  in  TEM  show  that  the  majority  of 
scattering  originates  from  a  constricted  volume,  called  the  effective 
scattering  volume.  This  is  true  for  a  region  of  homogeneous 
turbulence  and  is  anticipated  to  be  true  for  a  region  of 
inhomogeneous  turbulence,  although  no  calculations  have  been 
made.  Estimates  have  been  given  of  the  size  and  shape  of  the 
effective  scattering  volume  as  a  function  of  turbule  size  for  an 
experiment  conducted  in  homogeneous  turbulence  and  a  uniform 
atmosphere.  In  this  report,  homogeneous  turbulence  is  retained,  but 
an  upwardly  refracting  atmosphere  is  assumed,  resulting  in  a  shadow 
zone.  Inclusion  of  the  shadow  zone  boundary  further  limits  turbule 
sizes  and  locations  from  which  significant  signals  reach  the  detector. 
The  finding  is  that  large  turbules  are  less  effective  scatterers  and 
that  effective  scattering  volumes  are  large  enough  for  the  number  of 
small  turbules  to  be  large.  Implications  of  this  finding  are  discussed. 
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1.  Introduction 


Six  papers  were  presented  at  the  Battlefield  Atmospherics  Conferences  and 
other  conferences  that  deal  with  acoustic  scattering  by  atmospheric  turbu¬ 
lence.  The  sixth  paper  (Auvermann  and  Goedecke,  1995a)  contains  a  sum¬ 
mary  of  the  first  five  papers;  this  report  is  an  extension  of  that  paper  and  the 
fifth  paper  in  the  series  (Auvermann  et  al,  1994).  Here  we  refer  to  the  fifth 
paper  in  that  series  (Auvermann  et  al,  1994)  as  paper  A  and  the  sixth  paper 
(Auvermann  and  Goedecke,  1995a)  as  paper  B.  Experimental  evidence  shows 
that  atmospheric  turbulence  near  the  ground  is  neither  homogeneous  nor 
isotropic,  two  conditions  required  for  the  usual  statistical  model  of  turbu¬ 
lence  to  be  valid.  An  alternate  model,  termed  the  Turbule  Ensemble  Model 
(TEM),  consists  of  a  collection  of  turbules  of  different  sizes.  A  turbule  is  an 
isolated  inhomogeneity  of  either  temperature  or  velocity.  In  TEM,  the  scat¬ 
tering  pattern  of  individual  turbules  is  assumed  known.  The  analysis  pro¬ 
ceeds  by  the  assumption  of  a  distribution  function  for  the  sizes,  and  then 
location  of  the  turbules  of  each  size,  randomly,  within  the  atmospheric  re¬ 
gion  of  interest.  The  shadow  zone  signal  is  then  the  summation  of  the  con¬ 
tributions  from  each  turbule. 

Acoustical  signals  of  interest  to  the  Army  are,  in  general,  low  frequency. 
The  importance  of  this  is  that  wavelengths  are  large  compared  to  the 
dimensions  of  either  source  or  detector.  Therefore,  both  source  and  detector 
are  nearly  omnidirectional  and  thus  cannot  serve  to  define  a  scattering  vol¬ 
ume.  Paper  A  (Auvermann  et  al,  1994)  addressed  the  problem  of  determin¬ 
ing  the  volume  from  which  significant  scattering  can  occur  (called  the  effec¬ 
tive  scattering  volume  (ESV))  for  an  elementary  spherical  shell  volume. 
Paper  B  (Auvermann  and  Goedecke,  1995a)  addresses  the  determination  of 
ESV  for  a  homogeneous  turbulence  region  isolated  from  the  ground.  In  this 
report,  a  shadow  zone  boundary  (SZB)  is  introduced  into  the  scenario,  with 
the  remaining  properties  of  the  medium  assumed  to  be  homogeneous.  The 
full  power  of  TEM  for  addressing  scattering  from  anisotropic  inhomoge¬ 
neous  turbulence  has  not  been  exploited  as  yet.  The  signal  from  anisotropic 
inhomogeneous  turbulence  may  be  calculated  by  populating  the  ESV  with 
randomly  oriented  and  located  turbules  and  summing  the  scattered  power. 
Presently,  the  summation  process  is  carried  out  with  the  use  of  integrals 
over  mathematically  defined  distributions.  A  distinction  must  be  made  be¬ 
tween  isotropic  turbulence  and  isotropic  scattering.  Isotropic  turbulence,  in 
the  context  of  TEM,  refers  to  an  ensemble  of  turbules  that  has  uniformly 
distributed  orientation  vectors.  Such  an  ensemble  has  anisotropic  scatter¬ 
ing  properties  that  delimit  ESV  as  a  function  of  turbule  size.  Only  velocity 
turbulence  will  be  considered  in  this  report. 
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2.  Background 


In  this  section,  the  results  obtained  in  papers  A  and  B  are  recounted.  Paper  A 
was  the  first  step  in  the  process  of  using  the  scattering  properties  of  turbules 
to  define  ESV.  The  calculation  there  was  simply  of  the  total  scattering  cross 
section  of  the  turbules  in  a  spherical  shell  surrounding  the  detector.  The 
source  was  assumed  to  be  a  long  distance  away,  so  that  the  incident  signal 
was  a  plane  wave.  The  turbulence  had  parameters  appropriate  for  a  10-m 
height  above  the  ground.  However,  in  the  simplified  calculation,  the  source 
and  detector  were  assumed  to  be  far  enough  above  ground  as  to  be  in  a  free 
atmosphere.  Sound  speed  gradients  were  assumed  to  be  zero,  so  that  no 
shadow  zones  were  present.  The  largest  turbule  considered  had  a  10-m 
radius,  which  was  appropriate  for  the  10-m  height  of  the  scenario.  The 
frequency  and  wavelength  were  500  Hz  and  0.688  m,  respectively.  The  scat¬ 
tering  pattern  of  the  large  turbules  was  found  to  greatly  limit  the  region  of 
space  from  which  the  scattered  signal  could  reach  the  detector.  In  spite  of 
this  small  region  of  space,  the  scattering  properties  of  these  large  turbules 
dominated  the  scattering  cross  section.  This  scenario  (a  plane  wave  from 
the  source  in  a  medium  of  infinite  size)  cannot  be  extended  to  calculate  the 
signal  from  an  infinite  space,  because  the  scattered  signal  at  the  detector 
would  be  infinite.  Although  the  signal  from  more  distant  spherical  shells 
falls  off  as  1/r2,  the  total  volume  of  the  shell  increases  by  r2,  leaving  the 
signal  from  each  shell  the  same.  When  the  total  is  computed  for  an  increas¬ 
ing  distance  r,  the  total  increases  without  bound. 

In  paper  B,  the  distant  source  scenario  was  abandoned  in  favor  of  an  ex¬ 
perimental  scenario,  with  the  source  and  detector  finitely  separated.  The 
1/r2  loss  from  source  to  scatterer,  coupled  with  the  1/r2  loss  from  scatterer 
to  detector,  served  to  produce  a  finite  signal  at  the  detector.  ESV  was  found 
to  vary  from  4.4  m3  for  the  largest  turbule  to  2311  m3  for  a  turbule  of  0.44  m 
radius.  Figure  1  (taken  from  paper  B)  shows  the  intersection  locus  of  the 
scattering  volume  limit  with  the  vertical  plane  containing  the  source  and 
detector  for  three  different  turbule  sizes. 

Figure  2  depicts  the  scenario  of  this  report.  The  coordinates  are  x,  y,  and  z, 
with  the  origin  on  the  ground  beneath  the  source.  The  source  and  detector 
heights  are  hs  and  hD,  respectively.  The  wind  flows  from  the  direction  of  the 
detector  toward  the  source.  The  simplest  wind  speed  profile  will  be  used  to 
produce  a  shadow  zone,  namely,  a  uniform  gradient  model.  The  plan — the 
execution  of  which  is  discussed  in  section  5 — is  to  introduce  a  wind-driven 
shadow  zone  boundary  into  the  scenario  of  figure  2  and  then  calculate  the 
ESV  in  a  manner  similar  to  the  calculation  in  paper  B.  A  reasonable  wind 
speed  gradient  will  show  an  appreciably  sized  shadow  zone  over  the  ex¬ 
periment  distance  chosen  (320  m).  Two  matters  not  addressed  in  paper  B 
must  be  dealt  with  first:  (1)  calculating  parameters  (turbule  spacing  and 
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size  increment)  for  a  representative  turbulence  distribution  and  (2)  devel¬ 
oping  an  equation  for  the  SZB  surface.  These  matters  are  discussed  in  the 
next  two  sections. 

Symbols  for  some  of  the  variables  and  parameters  used  in  this  report  are 
summarized  in  the  following  list: 

a  =  turbule  characteristic  size  (m) 

c  =  sound  speed,  3.44  m  s-1 

Eg  =  total  irradiance  at  detector 

Cp  -  relative  irradiance  at  the  detector 

/  =  acoustic  wave  frequency  =  500  Hz 

6  =  spherical  coordinate  system  polar  angle 

(j)  =  spherical  coordinate  system  azimuthal  angle 


k 

l 

V 

X 

r 

Rs 

Rsd 

Rsp 

RpD 

CviV'l) 


vn{r,  6) 
V , 


X„ 

X,w 

X„ 


V 

Q 


2  n/k 

ratio  of  turbule  average  spacing  to  turbule  characteristic  size  n 

turbule  size  ratio  increment  parameter 

wavelength,  m  =  0.688  m 

total  number  of  turbules  of  the  largest  size 

spherical  coordinate  system  radial  distance 

source  irradiance  reference  distance  =  1.0  m 

source  detector  distance  =  320  m 

source  (point  P)  distance 

(point  P)  detector  distance 

velocity  turbule  total  scattering  cross  section 

turbule  flow  velocity  parameter 

turbule  velocity  function 

scattering  volume 

size  parameter  =  kn 

upper  limit  of  size  parameter  considered  =  91.325367 
maximum  in  the  size  parameter  integration 
minimum  in  the  size  parameter  integration 
scattering  angle 

magnitude  of  the  angular  velocity  vector 


3.  Turbule  Spacing  and  Size  Increment  Parameters 


T(n) 


The  turbule  spacing  and  size  increment  parameters  are  the  /  and  p,  respec¬ 
tively,  from  the  list  of  parameters  in  section  2.  The  size  increment  parameter 
will  be  determined  by  energy  conservation  considerations  for  a  turbule 
velocity  distribution  of  (Goedecke,  1992) 


v(f)  =  Qzxr  exp  (-  r2/a2)  , 
v„(r,  0)  =  (va0/a)  r  sin  (0)  exp  (-  r2/a2)  , 


where  the  rotation  axis  has  been  set  along  the  z-axis.  Because  of  rotational 
symmetry,  the  latter  expression  shows  the  dependence  on  the  spherical  co¬ 
ordinates  of  r  and  0.  The  total  kinetic  energy  of  the  turbule  T(a)  will  be  the 
integral  of  one-half  the  density  p  times  the  square  of  this  velocity: 


'pVaO 
v  2  ci2 


J  dr  r4exp  (-2  r2/a2)  J  d0[sin(0)]3J  d(j>  = 


n  pY  (5/2,  0,~) 


3x2 


1/2 


2  3 
VMa 


(2) 


where  Y  is  the  Generalized  Incomplete  Gamma  Function  (Wolfram,  1991). 
The  ratio  of  the  energy  of  the  largest  turbule  to  that  of  the  next  smaller  one 
using  the  former  velocity  relation  (Auvermann,  1994)  is  the  following: 


IK) 

T(„2) 


V 2  fl3 


2  3 
£>2  ^2 


11/3 


=  [exp  (p)] 11/3  -  2;  p  =  (3/11)  In  (2)  -  0.1890  .  (3) 


The  first  expression  in  equation  (3)  results  if  the  largest  turbules  break  up 
into  two  turbules  of  the  next  lower  size  and  split  energy  equally.  There  are 
49  size  classes  in  this  distribution,  from  the  largest  at  10  m  (denoted  by  the 
subscript  1)  to  the  smallest  at  1.146  mm  (denoted  by  the  subscript  49). 


Next,  recall  the  expression  for  the  velocity  structure  constant  (Goedecke  et 
al,  1998),  which  is 


jlu/3  (°'  °°)  =  dv  v 14/3  exP  (-  y2/ 2)  =  211/6r  (17/6, 0,  ~)  =  6.1455  . (4) 


For  fractal  scaling,  each  turbule  size  class  has  the  same  J  (Goedecke  et  al, 
1998).  N\/Vs  is  the  concentration  of  the  largest  turbules,  so  that  the  second 
bracket  is  l~3.  The  value  of  v-y  that  we  have  been  using  is  3.44  m  s-1.  To  obtain 
a  value  for  the  structure  parameter,  we  resort  to  the  literature  (Brown  and 
Clifford,  1976).  The  velocity  structure  constant  expression  is  given  as 

=  0.04  + 0.33  z~  2/3  ,  (5) 
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where  z  is  the  height  in  meters.  Substitution  of  10  for  z  gives  a  value  for  the 
velocity  structure  parameter  of  0.1111  m4/3-s“2.  Rewriting  equation  (4)  gives 
the  value  of  /: 

,  -  ( ft694%#--r  -  f  -<&69)(3'44)2(61455>)'/3  -  8.0148  .  (6) 

\  C;  fj  rt2/3  J  \ (0.1111)  (0.1890)  (10)2/3j 

If  the  centers  of  vortices  are  closer  than  2.8  times  the  diameter  of  an  indi¬ 
vidual  vortex,  the  mutual  distortion  can  be  expected  to  cause  the  vortices  to 
disintegrate  (Moore  and  Saffman,  1975,  eq  (13),  p  468).  This  number  of  2.8 
would  translate  into  a  value  of  5.6  if  the  spacing  were  compared  to  the  vor¬ 
tex  radius.  Thus,  the  value  of  8  from  equation  (6)  is  reasonable,  although 
the  exact  relation  of  the  vortex  diameter  from  Moore  and  Saffman's  work  to 
the  turbule  characteristic  size  (a)  of  the  model  is  not  known. 


4.  Derivation  of  Shadow  Zone  Boundary  Equation 


For  a  wind  speed  of  3.44  m  s-1  at  a  height  of  10  m  and  zero  at  the  ground, 
the  gradient  g(0)  is  0.344  s-1.  A  uniform  gradient  profile  declining  with  in¬ 
creasing  x  produces  a  circular  ray  path  (Pierce,  1989,  p  385),  whose  radius  is 
the  length  for  the  sound  speed  to  extrapolate  to  zero  or  1000  m.  Negligibly 
small  transverse  curvature  of  ray  paths  is  assumed,  so  that  a  ray  launched 
in  a  vertical  plane  making  an  angle  (3  with  the  x-z  plane  will  remain  in  that 
plane.  The  sound-speed  gradient  is  related  to  the  component  of  the  wind 
velocity  in  the  launch  plane  or  gQ 3)  =  g(0)  cos(/J).  The  SZB  is  the  limit  ray 
circle  that  contains  the  source  point  and  is  tangent  to  the  ground.  For  /?  =  0, 
the  equation  of  this  circle  is  given  in  the  following: 

(*S_rrf))  +{Z~Zco)  =  rc0'  rcO  =  C//&(0)'  Zc0  =  rc0^S~^s) 

Jrs  =  {rto-[rrt-(2-2<o)2]1/2>  •  <7> 

The  last  form  in  equation  (7)  is  the  useful  one,  since  the  height  of  the  SZB  at 
a  variable  z  will  later  determine  a  limit  for  an  integral.  For  the  circle  in  a 
plane  at  an  angle  /3  with  the  x-z  plane,  in  equation  (7),  replace  y(0)  with  g(fi) 
and  (z,  zc0)  by  analogous  distances,  say,  (ftp,  pcp).  Then,  replace  cos(j8)  and 
(ftp,  pcp)  with  the  appropriate  functions  of  the  variable  coordinates  (y,  z): 

{xS~rcp)2  +  [Pp-pcp)2  =  rcp;  Pp=(lJ2  +  z2)l/2 '  cos(P)  =  z/pp  . 

rcp=  c/(g(0)cos(/3)) ;  pcp=  (l  rcphs-h2sf/2  . 

xs  =  {rcp-[r2cp-{Pp-Pcp)2}V2} '  z>0  •  (8) 

=  0,  z  <0  . 

The  last  expression  for  xs,  where  z  is  less  than  or  equal  to  0,  recognizes  that 
propagation  in  the  negative  direction  (down  wind)  will  not  experience  a 
shadow  zone.  This  is  unimportant,  except  for  algorithm  purposes.  There  is 
also  an  SZB  associated  with  the  detector.  Signal  scattered  from  turbules  be¬ 
low  the  detector  SZB  will  not  reach  the  detector.  This  boundary  will  be  of 
the  same  form,  except  that  the  z  coordinate  will  turn  to  -zD,  and  then  the  zD 
will  be  replaced  by  (RSD  -  z) 

(■ XD-rca)2  +  {Pa-Pccc)2  =  r2ca '  Pa=\y2  +  {RSD-Z)2]l/2  '  COS  (a)  =  (RSD~z)/pa  . 

rca  =  c/(g(0)cos(a ));  pca=  (2 rcahD-h2D) 1/2  . 

xD={rca-{r2ca-{Pa-Pca)2}m} '  0<Z<RSD  (9) 

=  0,  z  >  Rsd  . 

The  integral  limit  mentioned  above  will  be  determined  for  a  y-z  plane  point 
(y,  z)  by  choice  of  the  larger  of  xs  or  xD.  The  larger  of  X5  or  xD  is  referred  to  as 
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5.  Scattered  Signal  at  Detector  Equation 


Since  we  make  the  assumption  that  the  only  influence  of  the  wind  gradient 
and  the  presence  of  the  ground  is  to  produce  the  SZB  of  section  4,  the  equa¬ 
tion  for  the  scattered  signal  at  the  detector  will  be  set  up  in  terms  of  the 
distances  and  angles  relative  to  the  source-detector  line  of  figure  2.  And 
since  we  also  assume  homogeneous,  isotropic  turbulence,  the  scattering  will 
have  rotational  symmetry  with  respect  to  this  line.  Imposition  of  the  SZB 
will  produce  only  a  rotational  angle  in  the  <j>  integration.  For  the  velocity 
distribution  of  equation  (22),  the  orientation  averaged  cross  section  is 
(Goedecke,  1992),  with  ca 0  as  the  asymptotic  sound  speed, 

<sv(¥,  X)  =  (f )  [sin(i^)cos(iy)]2exp{-22[l-cos(iy)]}  .  (10) 

The  total  signal  will  be  the  differential  signal  multiplied  by  a  differential 
volume  with  a  space  integration.  Including  integration  over  size  param¬ 
eter,  the  final  expression  is  in  equation  (8)  (Auvermann  et  al,  1994).  The 
differential  volume  is  a  partial  ring  around  the  source-detector  line. 


D  {  24  (ep- 1) 


nvx  a1  Rs\2  Cz 


)  fz  dZ  Jo  " dy  ^  ^  COS^  (Rsp  RpD )" 


[sin(^)cos(i/)]2  f  '"^/4/3exp{-22[l-cos(f)]}  , 

J  X„ 


(11) 


with  the  ring  ending  on  the  SZB  on  each  side.  By  symmetry  of  the  SZB  with 
respect  to  the  x-z  plane,  the  differential  volume  is  2  r,.(y,  z )  <p(y,  z)  dz  dy, 
where  rr  is  the  radius  of  the  ring  and  <j>  is  the  ring  angle  between  the  x-z 
plane  and  the  SZB  intersection  point.  In  equation  (11),  (r,.,  0,  R$p,  Rpd>  W)  are 
all  functions  of  (y,  z)  through  the  SZB  height  hz.  The  cos(y)  accounts  for  the 
inclination  of  the  source-detector  line  [tan(y)  =  (hD  -  hs)/Rso]-  Es  is  the  ref¬ 
erence  irradiance. 
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6.  Imposition  of  Significance  Criteria 


Expressions  for  the  relative  signals  of  interest  are 

^  i^i,' Zni,Yin,Xn,Xlu)  =  f  d  x  f  dz  f  dy\2  i'r  0](KSp  Rpp) 

JX “  ,  JZ"  J°  (12a) 

[sin(  i//)cos(  y/)]~  7 14/3  exp  {-  j2  1 1  -  cos(  y/)  j}  , 

1  (-oo,  oo,  oo,  o,  XJ  =  f  dz  f  %  [2  rr  0](/?s;,  R,,d)“2 

Joo  JO 

[sin(^)cos(^)|2  (*  ' '^J/ir14/'>exp{- j2|l-cos(v/)|}  ,  12b 


■f  -fJ(Z,,>Z,H>yH|/X„,X||t)l_ 

D  \  /(-<*>,  «,  0,  xj  j 


(12c) 


(Z„,  Z„„  Y„„  X,„  XJ  are  the  integration  limits  for  two  conditions.  The  first 
condition  is  that  for  all  space  for  which  these  limits  are  infinity,  zero,  or  X„. 
The  objective  of  this  report  is  to  find  limits  subject  to  the  condition  that  the 
signal  is  a  large  fraction  of  the  infinite  limit  signal,  equation  (12b).  Thus  the 
expressions  may  dispense  with  the  constant  in  front  of  the  integrals  in 
equation  (11).  The  relative  signal  cjj  is  defined  in  terms  of  the  integrals,  and 
the  significance  criterion  1  is  set  equal  to  it.  The  meaning  is  that  the  limits 
are  to  be  found  that  make  at  least  as  great  as  I,  which  is  thought  of  as 
being  a  number  such  as  0.99. 
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7.  Determination  of  ESV 


The  procedure  to  be  used  is  the  following.  The  reference  integral  (eq  12(b)) 
will  be  calculated  for  the  parameters  specified  in  this  report.  Then,  lower 
and  upper  limits  for  the  %  integration  (X„,  X,„)  will  be  determined  such  that 
the  integral  will  differ  in  each  case  by  no  more  that  0.002  from  the  reference 
integral,  with  all  other  limits  the  same.  Then  the  y  integral  will  be  evaluated 
as  a  function  of  (x,  2)  to  determine  a  table  of  Ym(x,  z)  values  such  that  this 
integral  differs  no  more  than  0.002  from  its  value  when  the  upper  limit  is 
infinity.  A  function  that  approximates  this  table  will  be  determined.  Finally, 
the  integral  with  these  limits  will  be  calculated  to  show  that  the  relative 
signal  is  greater  than  the  criterion  X.  These  criteria  are  summarized  in  table 
1.  Table  2  summarizes  the  calculated  results.  Table  3  has  been  constructed 
for  five  of  the  sizes  of  equation  (3).  From  table  2,  it  is  seen  that  the  conse¬ 
quence  of  this  rather  simple  criteria  scheme  was  that  Y,„  approximated  a 
cylinder,  and  the  total  X  was  too  large.  This  implies  that  a  more  sophisti¬ 
cated  scheme  in  which  the  x  dependence  of  the  z  limits  is  required  to  im¬ 
prove  the  accuracy  of  the  results.  The  SZB  is  likely  to  be  more  like  that  of 
figure  4  than  a  cylinder.  However,  important  conclusions  can  be  drawn 
using  the  results  of  table  3.  The  SZB  length  ls  has  been  estimated  from  other 
data,  and  the  concentration  Cx  has  been  calculated  for  an  l  (eq  (6))  of  8. 
Associated  with  the  table  are  some  numbers  calculated  for  the  second  en¬ 
try.  This  entry  is  the  maximum  contributor  to  the  total  signal.  Dividing  ls  by 
Nx  for  this  size  gives  an  average  spacing,  as  shown.  Dividing  the  average 


Table  1.  Selection  of 
significance  criteria. 


Goal:  X  =  0.99 

X  (X„)  =  0.998 
X  (Xm)  =  0.998 
X  (Z„)  =  0.998 
X  (Z,„)  =0.998 
X  (Ym)  =0-998 


Table  2.  Scattering 
volume  results. 


o,x„) 

=  0.456394 

■11 '  x,„) 

=  0.452400 

X 

=  0.9912 

X„ 

=  1.3249 

Xm 

=  52.4223 

Z„ 

=  44.5736  m 

z,» 

=  272.8049  m 

Y,„ 

=  816.4478  (x )' 
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Table  3.  Scattering 
volume  constituents. 


X 

ax 

r. 

h 

V, 

cx 

"x 

51.80 

5.67 

4.63 

40.00 

1.350x10" 

1.072  x  10-" 

1.440  x  lO"2 

16.67 

1.83 

20.47 

120.00 

7.898  x  104 

3.187  x  10-4 

2.517  xlO1 

7.83 

0.86 

55.12 

182.78 

8.723  x  10" 

3.071  x  10~" 

2.678  x  10" 

2.52 

0.28 

243.51 

220.49 

2.054  x  107 

8.897  x  lO'2 

1.827  x  106 

1.43 

0.16 

511.84 

228.23 

9.392  x  107 

4.768x10-' 

4.478  x  107 

X 

=  size  parameter 

Dominant  scatterer  has 

ax 

=  turbule  radius,  m 

size  parameter  =  16.67 

r. 

=  scattering  volume  radius,  m 

Spacing  in  ESV  =  4.767  m 

l 

=  scattering  volume  length,  m 

Wind  speed  =  3.44  m-s' 

V, 

=  scattering  volume,  m3 

Period  =  1.4  s 

C'z 

^x 

=  turbule  concentration,  nr3 
=  number  of  turbules 

Experimental  period  =  4  s 

spacing  by  the  wind  speed  gives  an  estimate  of  the  period  that  might  be 
expected  as  turbules  of  this  size  are  tracked  through  the  ESV.  This  period  of 
1.4  s  is  close  to  the  period  from  experimental  results.*  This  fact  is  important 
because  it  suggests  that  a  more  accurate  calculation  is  warranted.  However, 
a  more  accurate  calculation  will  need  a  realistic  estimate  of  the  size  param¬ 
eter  division  point  between  that  part  of  the  signal  that  can  be  obtained  by 
integration  (a  continuum  calculation)  and  those  parts  that  must  be  addressed 
as  discrete. 


*Shadow  zone  data  reported  by 
privately  by  the  authors  to  arrive  at 


Havelock,  Stenson,  and  Daigle  (1992)  were  examined 
this  estimate. 
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8.  Conclusions 


The  ultimate  utility  of  the  work  on  ESV  may  be  that  it  will  afford  a  physical 
explanation  of  the  shadow  zone  scattered  signal  fluctuations.  The  data  from 
table  3  of  section  7  lead  to  the  belief  that  we  can  conclude  that  this  variabil¬ 
ity  is  caused  by  the  motion  of  a  limited  number  of  moderately  sized  turbules 
through  the  ESV.  The  next  step  in  the  investigation  will  require  a  more  so¬ 
phisticated  method  for  determining  the  (z,  y)  limits  of  the  scattering  vol¬ 
ume.  Of  equal,  or  even  greater,  importance  will  be  the  determination  of  the 
minimum  population  of  the  scattering  volume  as  a  function  of  size  param¬ 
eter  that  will  result  in  insignificant  fluctuations.  This  minimum  population 
will  determine  the  turbule  size,  which  divides  the  sizes  that  must  be  summed 
over  from  the  sizes  that  can  be  integrated. 
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